Transitions in the wake of a flapping foil 
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We study experimentally the vortex streets produced by a flapping foil in a hydro dynamic tunnel, 
using 2D Particle Image Velocimetry (PIV). An analysis in terms of a flapping frequency-amplitude 
phase space allows to identify: 1) the transition from the well-known Benard-von Karman (BvK) 
wake to the reverse BvK vortex street that characterizes propulsive wakes, and 2) the symmetry 
breaking of this reverse BvK pattern giving rise to an asymmetric wake. We also show that the 
transition from a BvK wake to a reverse BvK wake precedes the actual drag-thrust transition and 
we discuss the significance of the present results in the analysis of flapping systems in nature. 



I. INTRODUCTION 

Control of vortices produced by flapping wings or fins 
to generate propulsive forces is the everyday task of birds, 
insects and swimming animals. Many studies of actual 
flapping extremities have been driven by the need for 
a better understanding of this form of propulsion with 
the ultimate goal of enhancing man-madepropulsive de- 
vices [see extensive reviews in U, d, d, 0, [5|. A wake 
vortex street with the sign of vorticity of each vortex 
reversed with respect to the typical Benard-von Kar- 
man (BvK) vortex street is an ubiquitous feature related 
to propulsive flapping motion. Such reverse BvK vor- 
tex streets have not only been observed in the wakes 
of swimming fish [see e.g. 0, 01 but also studied in de- 
tail through laboratory experiments with flapping foils 
& i^EEE^ [ill and numerical simulations 

0, H [li, [li, [23|. The mean velocity profile asso- 
ciated with these thrust-producing reverse BvK streets 
has the form of a jet and the stability properties of this 
average jet flow seem closely related to the propulsive 
efficiency of the flapping foil [21i] . 

Flapping systems, either natural or man-made, are of- 
ten discussed in terms of a single parameter, the Strouhal 
number, defined as the product of the flapping frequency 
/ and amplitude A divided by the cruising speed U, i.e. 
Sta = fA/U. Experiments with flapping foils [1, [22| 
have shown that propulsive efficiency peaks within a nar- 
row interval 0.2 < Sr^ < 0.4. Data compiled by [HI, [H] 
shows that observed cruise Strouhal numbers for a wide 
range of flying and swimming animals also lie on this in- 
terval, suggesting that, not surprisingly, natural selection 
has tuned animals for high propulsive efficiency. How- 
ever, less clear are the physical reasons that determine 
certain flapping configurations, with their associated vor- 
tex signature in the wake, to be the optimal choices for 
efficient thrust generation. 

Although three-dimensionality plays an important role 



in most flapping systems, a quasi-two-dimensional pic- 
ture of the wake often contains key dynamical elements, 
such as the creation and organization of vorticity, that 
are crucial in every system involving flapping foils as a 
means of producing propulsive forces [H, UM, HI, HE]- 
The goal of the present work is to explore these quasi-2D 
mechanisms using a simple but very well controlled ex- 
periment, in order to provide a definitive framework for 
the analysis of the transitions observed in the wake of 
flapping foils. 



II. EXPERIMENTAL SETUP AND 
PARAMETERS 

The setup consists of a pitching foil of aspect ratio 
4 placed in a hydrodynamic tunnel (see Fig. [1]). The 
experimental control parameters are the flow velocity in 
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FIG. 1: Schematic views of the foil profile and the hydro- 
dynamic tunnel. The foil chord c is 23mm and its span is 
100mm which covers the whole height of the 100 x 150mm 
section of the tunnel. The foil profile is symmetric, opening at 
the leading edge as a semicircle of diameter D = 5mm which 
is also the maximum foil width. The pitching axis is driven 
by a stepper motor. 
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the tunnel U, the foil oscillation frequency / and peak-to- 
peak amplitude A. The Reynolds number Re = UD/v^ 
where D is the foil width (or thickness) and v is the 
kinematic viscosity, was set to 255 in the experiments 
reported here. This corresponds to 1173 for the chord- 
based Reynolds number Rec = Uc/v^ which is an inter- 
mediate value in the 100 < Rec < 10000 range of 'inter- 
mediate Reynolds numbers' that characterize flapping- 
based propulsion in nature [5|. In addition, we deflne 
a dimensionless flapping amplitude and a Strouhal 
number Sr as 



Ad = A/D and Sr = fD/U . (1) 

Although natural vortex shedding exists for the steady 
flap at this Reynolds number even at zero-angle of in- 
cidence [35], no mode competition is observed in the 
strongly forced flapping regimes studied here. The flap- 
ping frequency used to deflne Sr is thus equivalent to the 
main vortex shedding frequency. In addition, we have 
chosen to use a flxed length scale (D) instead of the usual 
peak-to-peak amplitude (A) in the deflnition of Sr so that 
each degree of freedom of the flapping motion is repre- 
sented in one non-dimensional parameter. The {St^Ad) 
phase space thus deflned allows to give a clear picture of 
the different observed regimes. Measurements were per- 
formed using 2D Particle Image Velocimetry (PIV) on 
the horizontal mid-plane of the flap. PIV was performed 
using a LaVision system with an ImagerPro 1600 x 1200 
12-bit charge-coupled device (CCD) camera recording 
pairs of images at ~ 15Hz and a two rod Nd:YAG (15mJ) 
pulsed laser. Laser sheet width was about 1mm in the 
whole 100mm x 80mm imaging region. The time lapse 
between the two frames (dt) was set to 12ms. 



III. THE WAKE OF THE FLAPPING FOIL 

For a flxed Reynolds number, a point in the {St, Ad) 
parameter space deflnes thus the flapping conflguration 
and corresponds to a given vorticity signature in the 
wake of the foil. A parametric study as a function of 
the Strouhal number and the non-dimensional flapping 
amplitude allows to characterize the different regimes in 
the aforementioned (Sr, Ad) phase space and to pin down 
the transition zones. 



A. Symmetry properties 

In Fig. O we show snapshots of spanwise vorticity 
(left column) as well as average horizontal velocity flelds 
(right column), where the flapping amplitude Ad var- 
ied for a flxed value of the Strouhal number. The flrst row 
{Ad = 0.36) is the typical case of low-amplitude flapping 
which produces a forced wake with features resembling 
a natural Benard-von Karman (BvK) vortex street, but 



with shedding frequenc y lo cked-in to the flapping fre- 
quency [see also [23, m, [2§| . In particular, the vorticity 
engendered at the boundary layers on each side of the 
foil is rolled up in vortices that, in their downstream evo- 
lution, stay on the same side of the symmetry line of the 
wake, deflned by the axis of the foil when parallel to the 
free stream velocity. Accordingly, the mean flow in this 
case has a typical wake proflle with slower speed behind 
the obstacle. Increasing the flapping amplitude allows 
us to get a very clear view of the mechanism of rever- 
sal of vortex position: for amplitude Ad = 0.71 vortices 
of alternating signs align on the symmetry line of the 
wake (a 2S-type wake, using the nomenclature of (3Q|). 
whereas for = 1.07, the vortices shed on one side of 
the flap organize themselves in the wake on the other side 
of the symmetry line constituting the reverse BvK vortex 
street. The mean flow picture for this transition shows 
that the velocity deflcit in the wake proflle is almost com- 
pletely erased when vortices are aligned in the symmetry 
line {Ad = 0.71) prior to the apparition {dit Ad = 1.07 
and more clearly dit Ad = 1.77) of the jet proflle that is 
typically associated to the reverse BvK street [see also 
recent numerical simulations by 18]. 

Further increase in the flapping amplitude triggers a 
symmetry breaking in the reverse BvK wake (a feature 
previously observed for foils performing heaving motion 
[ill. fl3l. [31I . A strong dipolar structure that propagates 
obliquely to one side of the symmetry line of the wake 
is formed in each flapping cycle, while a much weaker 
single vortex is shed on the other side. In these asym- 
metric wakes, the quasi-two-dimensional nature of the 
flow is rapidly lost and the coherence of vortical struc- 
tures in the horizontal mid-plane where the observations 
are performed does not prevail for more than two cycles. 
However, the symmetry features of these wakes do seem 
to be deflned by a quasi-two-dimensional mechanism that 
depends on the initial conditions: the flrst dipole that is 
formed entrains fluid behind it, deflecting the mean flow 
in the wake and forcing the subsequent dipoles to fol- 
low the same path. This mechanism of dipole formation 
resulting in a deflected wake has also been observed in 
soap fllm experiments studying the wake of oscillating 
cylinders [32|, supporting the idea of a mainly 2D phe- 
nomenon. In the mean flow picture, the asymmetry of 
the wake can be clearly seen as a deflection of the peak 
in the jet proflle. 



B. The {St, Ad) phase space 

The experimental phase diagram is shown in Fig. [31 
where different symbols represent different types of wake. 
Experiments where a BvK-type wake was observed (□ 
symbols) occupy primarily the lower left region of the 
plot but extend to higher Strouhal numbers for the lowest 
amplitude tested. The transition from the BvK regime 
to the reverse BvK regime (blue line in Fig. [3]) can be 
univocally deflned in the {St, Ad) space when the posi- 
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FIG. 2: (Color) Instantaneous spanwise vorticity fields (left column, s ^) and mean flow (time averaged horizontal velocity, 
right column, m s~^) for fixed Strouhal and Reynolds numbers (Sr — 0.22 and Re = 255) and from top to bottom, for 
Ad = 0.36,0.71,1.07,1.77, and 2.8. The field of view (placed at mid-height of the foil) covers from -2D to 20D on the 
horizontal (streamwise direction x) and -8D to 8D in the vertical (crossstream direction y) where the origin is defined at the 
trailing edge of the flap at zero angle of attack. 
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FIG. 3: (Color) Ad vs. Sr map for Re = 255. Experimental 
points are labeled as □: BvK wake; ■: aligned vortices (2S 
wake); +: reverse BvK wake; A: deflected reverse BvK street 
resulting in an asymmetric wake. Blue line: transition be- 
tween BvK and reverse BvK. Green line: transition between 
reverse BvK and the asymmetric regime. The shaded area 
corresponds to the Sr^ = 0.3 ±0.1 interval, where Sr^ is the 
amplitude based Strouhal number SrA = St x Ad (see text). 



tion of vorticity maxima, at any given time and horizon- 
tal position, cross the symmetry axis of the wake. We 
remark that the transition Hne tends to an asymptotic 
value of Ad ~ 0.6 for Sr > 0.4 so that a threshold ampli- 
tude value exists for the production of a reverse street. 
The region where the reverse BvK regime is observed (+ 
symbols) is bounded on the other side by the transition 
to asymmetric regimes (A symbols) represented by the 
green line in Fig. [3l 



C. The drag— thrust transition 

Using the mean velocity fields to perform a standard 
momentum balance in a control volume enclosing the foil 
(see e.g. [33], pp. 349-351) we can obtain an indirect 
estimate of the mean drag. It is given by 



Fd = pUo j [Uo - u{y)]dy , (2) 

where p is the water density, Uo is the free stream ve- 
locity at the center of the tunnel upstream of the fiap 
and u{y) is a velocity profile measured in the wake. The 
position of this wake profile is set to X ^ 12 D down- 
stream from the trailing edge, where the wake can still 
be reasonably considered quasi-2D and across- wake pres- 
sure fiuctuations can be neglected. A drag coefficient 



surface Cd/Cdo obtained by interpolation of the mea- 
sured points in the (Sr, Ad) space are plotted as labeled 
iso-contours in Fig. [H The drag coefficient 



normalized by its value for the non-ffapping foil at zero 
angle of attack Cdoj allows to pin down the {St, Ad) 
curve where Cd changes sign, marking the transition be- 
tween drag and propulsive regimes (heavy black line in 
Fig. SI). A comparison of this zero-drag curve with the 
transition from a BvK vortex street to a reverse BvK 
pattern in the wake (gray line) shows that, for increas- 
ing ffapping amplitude or frequency, the reversal of the 
vortex street happens before the actual drag-thrust tran- 
sition in almost all the parameter range studied here. 
This shift shows that a region in the (Sr, A^) plane ex- 
ists where a reversed BvK pattern can be observed in the 
wake of the ffapping foil, but the relative thrust engen- 
dered by the ffapping motion is not yet enough to over- 
come the total mean drag. For Sr > 0.4, the two lines 
cross and the estimate of the drag-thrust transition stays 
actually slightly below the observed BvK-reverse BvK 
transition threshold. A possible explanation for this is 
that the quasi-2D hypothesis underlying the drag coeffi- 
cient calculation degrades faster the higher the Strouhal 
number, so that 3D effects might be cooperating to pro- 
duce this result that is unexpected from a 2D point of 
view. 

It should be noted that the approximations underlying 
the momentum balance used to obtain Eq. ([2]) make it a 
crude estimate for the mean drag on the ffapping foil. In 
particular, as soon as the ffow ceases to be parallel the 
momentum balance changes and Eq. ([2]) no longer holds. 
The mean drag in the regimes where the average wake 
is asymmetric and a mean lift force is likely to be non- 
negligible, will thus be poorly estimated by Eq. ([2]). We 
nonetheless verify a posteriori that the expected drag- 
thrust transition does not overlap the region of the pa- 
rameter space where asymmetric wakes develop, so that 
the prediction should be reasonable. The sensitivity of 
the drag coefficient calculation to the x-position of the 
wake proffle, which deffnes the downstream boundary of 
the control volume, was tested performing the momen- 
tum balance at different horizontal stations in the wake. 
Thus, the gray 'error band' around the Cd = transition 
line in ffgure [H represents the ffuctuation of the transi- 
tion line when the calculation of is performed using 
a wake proffle measured at X ± AD (i.e. representing an 
overall span in the x direction of ^ 40% of the length of 
the PIV measurement window) jH]. 

IV. CONCLUSIONS 

The experimental study of the wake vortices produced 
by a high-aspect-ratio pitching foil presented in this pa- 
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FIG. 4: Contours of a mean drag coefficient Cd/Cdo surface 
estimated using a momentum-balance approach (see text). 
The black line corresponds to Cd = where the estimated 
drag-thrust transition occurs. The shaded area represents the 
estimated error for the Cd = curve due to sensitivity on the 
choice of the control volume. The gray line is the transition 
from BvK to reverse BvK (blue line in Fig. [3]). The dashed 
line corresponds to Sr^ = 0.3. 

per gives evidence characterizing two key dynamical fea- 
tures relevant to wake vortex systems engendered by flap- 
ping motion: flrst, the transition from the well-known 
Benard-von Karman (BvK) wake to the reversed vor- 
tex street that signals propulsive wakes, and second, the 
symmetry breaking of this reverse BvK pattern giving 
rise to an asymmetric wake. Albeit the inherently-three- 



dimensional nature of the wake dynamics resulting from 
the flnite span of the foil in the hydrodynamic tunnel, the 
quasi-two-dimensional picture studied here by focusing 
on the near wake captures the essential of the vorticity 
dynamics because the main vortex shedding is induced 
when the foil motion changes direction [see also S, 
where 2D theory and simulations show good agreement 
with experimental studies on flapping wings]. 

In conclusion, we can come back to the flapping-based 
propulsive regimes encountered in nature, where flsh-like 
swimming and flapping flight occur for Sr^ = 0.3 ± 0.1. 
In terms of the parameters used in this paper, Sr^ = 
St X Ad J that the previous interval represents a re- 
gion bounded by hyperboles in the (Sr, Ad) phase space 
(shaded area in Fig. [3]). Our flrst remark is that the drag- 
thrust transition shown in Fig. [4] compares well with the 
Sta = 0.3 curve (dashed line) that characterizes animal 
propulsion. Within the limits imposed by the quasi-2D 
picture analyzed in the present work, this suggests that 
the conclusions obtained here should be relevant as in- 
terpretation tools in the study of natural systems. In 
addition. Fig. [3] also shows that the region deflned by 
Sr^ = 0.3 ±0.1 overlaps not only with the reverse BvK 
regime, but also with the asymmetric regime. We may 
conjecture that animals using flapping-based propulsion 
must either exploit the creation of asymmetric wakes as 
part of their manoeuvering techniques or, when cruising, 
avoid flapping regimes where the symmetry breaking of 
the reverse BvK street will occur. 
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